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ABSTRACT 

We quantize the Chern-Simons-Proca theory in three dimensions by using the 
Batahn-Tyutin Hamiltonian method, which systematically embeds second class con- 
straint system into first class by introducing new fields in the extended phase space. As 
results, we obtain simultaneously the Stiickelberg scalar term, which is needed to can- 
cel the gauge anomaly due to the mass term, and the new type of Wess-Zumino action, 
which is irrelevant to the gauge symmetry. We also investigate the infrared property 
of the Chern-Simons-Proca theory by using the Batalin-Tyutin formalism comparing 
with the symplectic formalism. As a result, we observe that the resulting theory is 
precisely the gauge invariant Chern-Simons-Proca quantum mechanical version of this 
theory. 
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1 Introduction 



The Dirac method has been widely used in the Hamiltonian formahsm [1] to quan- 
tize second class constraint systems. However, since the resulting Dirac brackets are 
generally field- dependent and nonlocal, and have a serious ordering problem between 
field operators, these are under unfavorable circumstances in finding canonically con- 
jugate pairs. On the other hand, the quantization of first class constraint systems 
[2,3] has been well appreciated in a gauge invariant manner preserving Bccci-Rouct- 
Stora-Tyutin (BRST) symmetry [4,5]. This formalism has been extensively studied by 
Batalin, Fradkin, and Tyutin [6,7] in canonical formalism, and applied to various mod- 
els [8-10] obtaining the Wess-Zumino (WZ) action [11,12]. Recently, Banerjee [13] has 
apphed the Batahn-Tyutin (BT) Hamiltonian method [7] to the abehan Chern-Simons 
(CS) field theory [14-16]. As a result, he has obtained the new type of an abelian WZ 
action, which cannot be obtained in the usual path-integral framework. Very recently, 
we have quantized the nonabelian case by generalizing this BT formalism [17]. As 
shown in these works, the nature of second class constraint algebra in the original the- 
ories originates from the symplectic structure of CS term, not due to the local gauge 
symmetry breaking. Banerjee, Ghosh, and Banerjee [18] have also considered a massive 
Maxwell theory. As a result, the extra field in this approach has identified with the 
Stiickelberg scalar. We have also quantized the abelian self-dual massive theory by 
using this formalism, which interestingly produces both the Stiickelberg scalar and the 
new type of WZ [19]. There are some other interesting examples in this approach [20]. 

On the other hand, three-dimensional Chern-Simons gauge theories have been at- 
tracting much attention because these play an important role in the present devel- 
opment of the quantum Hall effect [21] and the string theory [22]. The quantum 
mechanical version of the CS field theory has been studied by Jackiw and collabo- 
rators through the phase space reductive limiting procedure [23]. Recently, Baxter 
[24] has described a simple (2-f l)-dimensional system which allow in principle as the 
experimental verification of the CS feature by introducing the Rontgen energy term 
[25]. 

In the present paper, we shall apply the BT Hamiltonian method [7] to the Chern- 
Simons-Proca (CSP) theory [26] revealing both the Stiickelberg effect [27] and the CS 
effect [13,17]. In Sec. 2, we apply the BT formalism to the CSP theory in three 
dimensions which is gauge non-invariant. By identifying the new fields p and A with 
the Stiickelberg scalar and the WZ scalar, respectively, we obtain simultaneously the 
Stiickelberg scalar term related to the explicit gauge-symmetry-breaking mass term 
and the new type of WZ action related to the symplectic structure of the CS term. 
In Sec. 3, we also investigate the quantum mechanical version of the CSP theory by 
using the BT formalism comparing with the symplectic formalism [29], which is the 
improved version of the Dirac method, and, in particular, very effective for the first- 
order Lagrangian. As a result, we observe that the resulting theory is just the gauge 
invariant CSP quantum mechanical model. 
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2 The Chern-Simons-Proca Theory 



We consider the abelian CSP model [26] 

1 



1 



(1) 



by using the BT formahsm. Note that this action has an exphcit mass term, which 
breaks the gauge symmetry as the case of the Proca model [18], and also the CS term, 
which has a different origin of the second class constraint system. Consequently, this 
action represents the second class constraint system combined with two effects, which 
can be easily confirmed by the standard Dirac analysis [1]. There are three primary 
constraints. 



^0, 

-KtijA^ ~ 0, 



1,2), 



and one secondary constraint, 

CU3 = rri^A^ - Keijd'A^ ^ 0, 
which is obtained by conserving fig with the total Hamiltonian, 

where is the canonical Hamiltonian, 



KSijA^'d'A^ + -mmA'f - (^0)2} 



(2) 



(3) 



(4) 



(5) 



and u^, are Lagrange multipliers. No further constraints are generated via this 
procedure. We find that all constraints are fully second class. In order to carry out 
the simple algebraic manipulation, it is, however, essential to redefine 013 by using 
as follows 



^3 = ^3 + d'VLi 



(6) 



although the redefined constraints J1q(q; = 0, 1,2,3) are still completely second class. 
Otherwise, one will have a complicated constraint algebra including the derivative 
terms, which is difficult to handle. Then, the modified constraint algebra is given by 



{Qa{x)Mp{y)} 
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which reveals the simple second class nature of the constraints Q^- 

In order to convert this system into first class, the first objective is to transform 
Qa into the first class by extending a phase space. Following the BT approach [7], 
we introduce new auxiliary fields and assume that the Poisson algebra of the new 
fields is given by 

{^'^{x),^^{y)} = u;'^f'{x,y), (8) 

where o;"^ is an antisymmetric matrix. Then, the modified constraint in the extended 
phase space is given by 



(9) 



n=l 



satisfying the boundary condition, 0^(71^, A'^, 0) — fla- The first order correction term 
in the infinite series [7] is given by 

= j <fyX^p{x,y)^P{y), (10) 

and the first class constraint algebra of Vta requires the condition as follows 

^apix, y) + j ^'^ X^^{x, w)u'^\w, z)Xsp{z, y) = 0. (11) 

Among the solutions satisfying with the conditions (8) and (11), we take a simple 
solution as follows 



u"^{x,y) 



Xapix.y) 



/ 1 \ 

0-10 

10 

V -1 / 

/ m \ 





V m y 



5'^{x-y), 



5\x-y). 



(12) 



There is an arbitrariness in choosing o;"^, which would naturally be manifested in Eq. 
(12). This just corresponds to the canonical transformations in the extended phase 
space. However, as has been shown in other calculations [13,17], this choice of Eq. (12) 
gives the remarkable algebraic simplification. 

Using Eqs. (9) and (12), the new set of constraints is found to be 



which are strongly involutive. 



VLq + m<I>°, 

= 0. 



1,2), 



(13) 
(14) 
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As a result, we have fully first class constraints in the extended phase space by applying 
the BT formalism systematically. Observe further that only the Q^^^ contributes in the 
series (9) defining the first class constraint. All higher order terms given by Eq. (9) 
vanish as a consequence of the choice Eq. (12). Recall the are the new variables 
satisfying the algebra (8) with c<;"^ given by Eq. (12). 

The next step is to obtain the involutive Hamiltonian, which naturally generates 
the secondary Gauss constraint in the extended phase space. It is given by the infinite 
series [7], 



(15) 



n=l 



satisfying the initial condition, H{7rn, = He- The general solution [7] for the 

involution of H is given by 

iy(n) ^ _i / d^xcfyd'z ^''(x)u;^p(x, y)X^'^(y, z)G^;'-'\z), (n > 1), (16) 

where the generating functionals G^"^ are given by 

= {Q^^\H(-% + {n(^\H^--'^o (n>l), (17) 

where the symbol O in Eq. (17) represents that the Poisson brackets are calculated 
among the original variables, i.e., O = {7ifj,,A^). Here, w^/? and X°'^ are the inverse 
matrices of co"^ and X^p respectively. Explicit calculations of G^^ yield. 



= m%A\ 

which are substituted in Eq. (16) to obtain H^^\ 



m^^diA' + ^eij$M^ + ^f^^'diA^ - $^(mA° - — e^aM^) 



(18) 



(19) 



This is inserted back in Eq. (17) to deduce G^^^ as follows 



(1) 



(1) 



G 



2 



m 



m 



which then yield H^'^^ from Eq. (16), 



1 777 777 ^ 



2k 



1 



m 



(20) 



(21) 
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Since G^^^ = (n > 2), the final expression for the involutive Hamiltonian after the 
n — 2 finite truncations is given by 

H = + H^^^ + H^^\ (22) 

which is strongly involutive, 

{Q^,H}^0. (23) 

According to the usual BT formahsm, this formally completes the operatorial conver- 
sion of the original second class system with Hamiltonian He and constraints Qa into 
the first class with Hamiltonian H and constraints Qa- 

However, before performing the momentum integrations to obtain the partition 
function in the configuration space, it seems appropriate to comment on the strongly 
involutive Hamiltonian. If we directly use this Hamiltonian, we can not naturally 
generate the first class Gauss' law constraint from the time evolution of the primary 
constraint Qq, which is the first class. Therefore, in order to avoid this problem, we use 
the equivalent first class Hamiltonian without any loss of generality, which only differs 
from the involutive Hamiltonian (22) by adding a term proportional to the first class 
constraint Cl^ as follows 

H'^H+—<^^Q3. (24) 

m 

Then, this desired Hamiltonian H' automatically generates the Gauss' law constraint 
such that {fio, H } = Q^. Note that when we act this modified Hamiltonian on physical 
states, the difference with H is trivial because such states are annihilated by the first 
class constraint. Similarly, the equations of motion for observable {i.e. gauge invariant 
variables) will also be unaffected by this difference since ^3 can be regarded as the 
generator of the gauge transformations. 

We now derive the Lagrangian, which will include both the Stiickelberg effect and 
the CS effect, corresponding to the Hamiltonian (24). The first step is to identify 
the new variables $° as canonically conjugate pairs in the Hamiltonian formalism as 
follows ^ ^ 

= (mp, -^7Tx, —TTp) (25) 

K TTl 

satisfying Eqs. (8) and (12). The starting phase space partition function is then given 
by the Faddeev formula [28] , 

r ^ 

where 

S' ^ J d^x (7r^> + ttaA + TTpp - H') (27) 

with the Hamiltonian density li! corresponding to i?', which is now expressed in terms 
of {p, TTp, A, tta} instead of The gauge fixing conditions Va may be assumed to be 
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independent of the momenta so that these are considered as the Faddeev-Popov type 
gauge conditions [28] . 

Next, we perform the momentum integrations to obtain the configuration space 
partition function. The ttq, tti, and 712 integrations are triviaUy performed by exploiting 
the deha functions S{flo) = ^{^o + m^p), S{fli) = 5{7ri + + tta), and 
5(^2) = ^{7^2 — f + i^X), respectively. After exponentiating the remaining delta 
function 5(^3) = S{—Keijd^A^ + diirx + Kd2\ + m^v4° + tt^) with Fourier variable ^ as 
5(^3) = /P^e-^/'^'^^^^ and transforming A^ + we obtain the action as follows 



S = j d^x {-]^Ke^,pA''d''AP + ]:m^A^'A, 



1 77?"^ 

+ pi-m" (AO + ^) - m^diA' - -m^did'p + m^diX Sstta] 

2 K 

+ 7r,[A -A'- "^A' + d'A' - In,] - ^m^f }, (28) 
where the over dot means the time derivative, and the corresponding measure is given 

by 

3 

[Dp] = VA'^VXVTrxVpVTTpV^ J] ^(r/3[^° + A\ A, p])det \ {Q^, F^^} | , (29) 

13=0 

where A^ — > ^4°+^ transformation is naturally understood in the gauge fixing condition 

r/3. 

Note that the original theory is easily reproduced in one line, i.e., if we choose the 
unitary gauge 

Fa = {p,7T\, X,7lp), (30) 

and integrate over ^. Then, one can easily realize that the new fields are nothing but 
the gauge degrees of freedom, which can be removed by utilizing the gauge symmetry. 

Now, we perform the Gaussian integration over TTp. Then all terms including ^ 
in the action are canceled out, and integrating over ttx the resultant action is finally 
obtained as follows 

S — Sst + Snwz + Sb 

Sst = Jd'x{-^Ke,,pA>^d^Af + ^m'{A^ + d,pf}, 
2 2 
Snwz = J ^'^^{^[A + Fqi + ^(^ + ^sp)]' 

+ nX[Fo2-—{A, + ^^p)-'^X]}, 
Sb ^ - J (fx d^im^pA^") (31) 
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where F^i, — dfj,Ai, — di,Aij,. Note that Sst is an expected Stiickelberg scalar term, which 
is needed to cancel the gauge anomaly due to the explicit gauge-symmetry-breaking 
mass term [18], Snwz is the new type of WZ term due to the symplectic structure of 
the CS term [13,19], which is irrelevant to the gauge symmetry, and Sb is the boundary 
term, which is also needed to make the second class system into the first class. The 
corresponding Liouville measure just comprises the configuration space variables as 
follows ^ 

[Dfj,] = VA'^VXDpVi n ^(r/3[A° + e, A\ A, p])det \ {h^, Vp} \ . (32) 

/3=0 

This action S is invariant up to the total divergence under the gauge transformations 
as 6A^ = d^A, 6p = —A, and S\ = 0. 

Note that starting from the action (31) with the boundary term Sb, we can easily 
reproduce the same set of all first class constraints Cla, and the Hamiltonian such that 

2 

/7TI 
(fx [nxdiAo TTA^a + rrfpdiA' + KXd2Ao + rrfXAi 

+ - ^n.d^p - lm%pd^p + m'Xd.p + \m'\' + ^^ttJ. (33) 

Then, if we add a term proportional to the constraint ^^3, i.e., ^TTpl^s, which is trivial 
when acting on the physical Hilbert space, to the above Hamiltonian (33), we can obtain 
the original involutive Hamiltonian (22). Furthermore, this difference is also trivial 
in the construction of the functional integral because the constraint ^3 is strongly 
implemented by the delta function 5{VL^) in Eq. (26). Therefore, we have shown 
that the constraints and the Hamiltonian following from the action (31) arc effectively 
equivalent to the original Hamiltonian embedding structure. As results, through the 
BT quantization procedure, we have found that the Stiickelberg scalar p is naturally 
introduced in the mass term, and this p as well as the WZ scalar A is also included in 
the new type of WZ action. 

We also note that if we ignore the boundary term Sb in this action, we cannot 
directly obtain the involutive first class Hamiltonian as the case of the Proca theory 
explained in Ref. [18] because this boundary term plays the important role in this 
procedure. 

Finally, note that in the trivial limit k ^ m, the action (31) exactly reduces to 
the self-dual massive theory having all the first class constraints, which has recently 
been derived in Ref. [19]. The hmit m — > is non-trivial because the action has 
the m"^ term. In fact, we can easily find that the auxiliary field tta is not well- 
defined in the action for the case of this limit because the ttx contains the term, 
i.e., Tlx = "%(^ + -foi) + ^(^2 + d2p). Therefore, we have to pay attention to when 
the momentum integrations are performed in the below of Eq. (28). Avoiding this 
situation, if we simultaneously take the limit (A-|- Fqi) —>■ with the m ^ 0, we obtain 
the delta function 5{X + Fqi) in the measure part when we perform the momentum 
integration over tta resulting to the case of the pure CS action. Then, one can finally 
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find that the CSP theory exactly reduces to the pure CS case having all the first class 
constraints [19]. 



3 The Chern-Simons-Proca Quantum Mechanics 

3.1 The Symplectic Quantization of the CSP Quantum 

Mechanics 

Let us briefly discuss the symplectic quantization, which is very effective for the 
first-order Lagrangian [29], of the CSP quantum mechanics. We start the following 
infrared limit action [23], which is already first-order, of the CSP theory with the 
Coulomb gauge 

dt -Keijq'q^ - -mq'q' , (34) 



where ei2 = e^^ = 1. Following the symplectic formalism, we first rewrite the action as 
follows 

(35) 



^ldt[-^q^' + ^q'q'-H('^ 



where H^'^^ = yH^^Y + (9^)^] is the usual canonical Hamiltonian and the superscript 
denotes the number of iterations [29]. Note that since the action (35) is already first- 
ordered from the start, we do not need to introduce auxiliary fields such as conjugate 
momenta. Then, we set symplectic variables = (?^,?^), and symplectically conju- 
gated momenta a^"'' = (— f q'^, f Q'^)- 

Now, symplectic 2- form matrix fij, which consists of the essential part for finding 
the generalized brackets of the symplectic formalism, is obtained as follows 



■lij Q^{0)i Q^{0)j I -1 i ^ ^ 



having the antisymmetric property. This matrix is not singular, and thus has the 
inverse as follows 



/<»'•' = 4 ( -1 ) ' 



(37) 



which gives the generalized symplectic brackets 



{q\q'}^--e''. (38) 

These are exactly same as the Dirac brackets [23] when we analyze the system through 
the usual Dirac's method. 

It is appropriate to comment on the symplectic formalism [29]. In general, the sym- 
plectic 2-form matrix is not invertible at the first stage of iterations. Then, we can find 
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some zero modes, which are related to generate the constraints in the symplectic for- 
mulation, and incorporate them into the canonical sector with some auxiliary variables 

to find the nonvanishing symplectic 2-form matrix. If we can find the invertible sym- 
plectic 2-form matrix at the finite stage of iterations, the inverse of the matrix gives the 
generalized brackets, which are equivalent to the usual Dirac brackets. However, when 
we can not find the invertible matrix even at the infinite stages of iterations, we can 
say the system has a gauge symmetry and use the zero modes to obtain the concrete 
rules of transformations [29]. Especially in the case of the CSP quantum mechanical 
model, wc have the symplectic 2-form matrix at the first stage of iterations. Thus the 
system has no constraints in the symplectic quantization formalism, while this system 
has second class constraints [23] in the standard Dirac formahsm as follows 

= + ^i^eijq^ ^0 (i = 1, 2), (39) 

where Pi — = —\^ijcf are canonical momenta. 

Now, using the generalized brackets, we can easily find the Hamilton equations as 
follows 

= = --eV- (40) 

K 

These equations give a single simple harmonic oscillator 

7(^) = g^ + ^g2 = ^a;7(^), (41) 

with the frequency a; = — as usual. In the next section, we will show that starting the 
gauge non- invariant action (34), we can obtain the gauge invariant version describing 
the simple harmonic oscillator by the BT formalism. 



3.2 The BT Quantization of the CSP Quantum Mechanics 

Now let us analyze the action (34) in the BT quantization as in the previous Section 
2. The first observed fact through the usual Dirac's procedure [1] is that the action 
represents a second class constraint system, i.e., there are two primary constraints (39), 
and no further constraints are generated through the time evolution of these constraints 
with the total Hamiltonian, 

Ht^H, + u%, (42) 
where is the canonical Hamiltonian, 

777, 

H, = p,ct-L, = -q'q\ (43) 
and are Lagrange multipliers. Then, the constraint algebra is given by 

(44, 
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which reveals the simple second class nature of the constraints Qj. 

In order to convert this system into first class, the first objective is to transform 
into the first class by extending the phase space. Following the BT approach [7], we 
introduce the matrix (8) as follows 

= ( J "o' ) • («) 

Then the other matrix X^j in Eq. (10) is obtained by solving Eq. (11) with the Aj^- 
given by Eq. (45), 

There is an arbitrariness in choosing u'^^ , which would naturally be manifested in Eq. 
(45) as explained in the Section 2. 

Using Eqs. (9), (45) and (46), the new set of constraints is found to be 

= + (47) 

which are strongly involutive, 

{a„n,} = o. (48) 

The next step is to obtain the involutive Hamiltonian. The generating functions 
G\^^ are obtained from Eq. (17). It is noteworthy that there are only two terms 
and in the expansion (47) due to the choice (45) and (46). Explicit calculations 
yield, 

= -mq\ (49) 
which are substituted in Eq. (16) to obtain H^^\ 

H^'^ = ^ei.^W- (50) 

This is inserted back in Eq. (17) to deduce Cf'^ as follows 

G« = ^6,$^ (51) 

which then yield H^^^ from Eq. (16), 

H(^) = H^i^i, (52) 

Since ^ = (n > 2), the final expression for the involutive Hamiltonian after the 
n — 2 finite truncations is given by 

H^Hc + H^^^ + = ^[q'q' + -^Sij^'q^ + -^W]. (53) 

2 y/K K 
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which, by construction, is strongly involutive, 

{ni,H} = 0. (54) 

This completes the opcratorial conversion of the original second class system with 
Hamiltonian He and constraints Qi into the first class with Hamiltonian H and con- 
straints fli- 

We now derive the gauge invariant Lagrangian corresponding to the Hamiltonian 
(53). The first step is to identify the new variables as canonically conjugate pairs 
in the Hamiltonian formalism, 

^- = {^P^,^X), (55) 

satisfying Eqs. (8), (45) and (46). The starting phase space partition function is then 
given by the Faddeev formula, 

Z= f Vq'VpiVPxVX ll6{ni)6{Tj)det \ {ni,Tj} \ e'^ , (56) 

id 

where 

jdt {p^q' + PxX - h) . (57) 

As similar to the Section 2, the gauge fixing conditions Fj may be assumed to be 
independent of the momenta so that these are considered as the Faddeev-Popov type 
gauge conditions. 

Next, we perform the momentum integrations to obtain the configuration space 
partition function. The pi, and p2 integrations are trivially performed by exploiting 
the delta functions 5{Vti) = 5(pi + fg^ + Px), and ^(^^2) = 6{p2 - ^q^ + nX), 
respectively. Then, we obtain the action as follows 



f 1 . . 1 . . 
S ^ J dt {-KSijq'q^ - -mq'q' 



1 . rji jji 

+ X{-Kq^ + mq'--mX) + Px{-q'+X--q''-—Px)}, (58) 
and the corresponding measure is given by 

[P/x] = Vq'VXVPx n S{Tj)det \ {Qi, Tj} \ . (59) 

Note that the original quantum mechanical model is also reproduced when we choose 
the unitary gauge such that Fq, = {X, Px) as the case of the CSP theory. 

Now, we perform the Gaussian integration over Px- Then, the resultant action is 
finally obtained as follows 

S — So + Snwz ; 

^2 



^NWZ 
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We can rewrite the above action with the boundary term more compactly as follows 

s = jdt{^{e-xr-^{q'-xr}+ss ■, 

Sb = + (61) 

The corresponding Liouville measure just comprises the configuration space variables 
as follows 

[Dfj] = Vq'VX II 6{Tj)det \ {tli, Tj} \ . (62) 
hj 

The action (61) is invariant up to the total divergence under the transformation 
Sq^ — e(t) and SX — e(t) , which are just the gauge transformations of the CSP quantum 
mechanical model. Note that the variable has not appeared in Eq. (61) except 
the boundary term Sb- Furthermore, the action except the term Sb is just a usual 
harmonic oscillator having the frequency a; = ^ as in the previous section when we 
define a quantity such that j{t) = — X. Since S^{t) is invariant under the above 
transformations, it is a physical quantity. As a result, starting from the gauge non- 
invariant system (34), we obtain the gauge invariant version describing the harmonic 
oscillator in the BT formalism. 

Finally, we would like to comment that the gauge invariant action (61) is not 
separated into the original action So and the new type of the WZ action Snwz- This 
is because X is nothing but the gauge degree of freedom. 

4 Conclusion 

In conclusion, we have applied the Batalin-Tyutin method, which converts the 
second class system into the first class, to the CSP theory and the quantum mechan- 
ical version of this theory. For the CSP case, we have shown that if we ignore the 
boundary term in action (31), the direct connection with the Lagrangian embedding 
of Stiickelberg scalar can be made by explicitly evaluating the momentum integrals 
in the extended phase space partition function using the Faddeev-Popov-like gauges, 
and identifying the extra field p introduced in our Hamiltonian formalism with the 
conventional Stiickelberg scalar needed to cancel the gauge anomaly due to the mass 
term. We have also obtained a new type of WZ action Snwz containing the WZ scalar 
A, which is irrelevant to the gauge anomaly. Furthermore, we should also keep the 
boundary term Sb- Otherwise, we cannot reproduce the original first class system. 
Note that the Stiickelberg scalar p is also included in S^wz in order to maintain the 
gauge invariance of the Snwz related to the CS effect in the action (31). 

On the other hand, we have observed that the infrared hmit of the CSP theory 
is precisely the gauge invariant CSP quantum mechanical model by using the BT 
formalism comparing with the symplectic formalism. Even though we can find a har- 
monic oscillator solution by solving the Hamilton equations of motion after applying 
the standard Dirac method, we have applied the symplectic method since it is more 
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intuitive when we directly find the generahzed brackets from the Hamilton equation 
fij^^ — ^H{^) for the first-order system. In other words, if fij has an inverse, the 

Hamilton equation is easily obtained through = {^*, H{^)} — {^', ^^} ^g^P , where 
{Cj ^"'} is the generalized bracket. Furthermore, by applying the BT formalism, we 
have firstly realized the harmonic oscillator at the action level, which is manifestly 
gauge invariant in the CSP quantum mechanical model. 
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